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We investigate the processes leading to phase-space equilibration of parton distributions in nu- 
clear interactions at collider energies. We derive a set of rate equations describing the chemical 
equilibration of gluons and quarks including medium effects on the relevant QCD transport coeffi- 
cients, and discuss their consequences for parton equilibration in heavy ion collisions. 
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I. INTRODUCTION 



The mechanisms that precede the formation of a thermalized, locally deconfined plasma of quarks and gluons in 
relativistic nuclear collisions have recently attracted considerable interest, because it was noticed that the preequilib- 
rium phase may influence the yield of certain quark-gluon plasma signals, such as lepton pairs and hadrons containing 
heavy quarks. The space-time evolution of quark and gluon distributions has been investigated in the framework of the 
parton cascade model jij. This model is based on the concept of the inside-outside cascade ^-[| and evolves parton 
distributions by Monte-Carlo simulation of a relativistic transport equation involving lowest-order perturbative QCD 
scattering and parton fragmentation. Numerical studies J5|,|| have shown that phase-space equilibration of partons 
occurs over a period of 1-2 fm/c, and is initially dominated by gluon-induced processes. 

From these investigations, and from more schematic considerations, a picture involving three distinct stages of 
parton evolution has emerged : 

(1) Gluons "thermalize" very rapidly, reaching approximately isotropic momentum space distributions after a time 
of the order of 0.3 fm/c. 

(2) Full equilibration of the gluon phase space density takes considerably longer. 

(3) The evolution of quark distributions lags behind that of the gluons, because the relevant QCD cross sections 
are suppressed by a factor 2-3. 

Although this picture emerges from the numerical simulations of the parton cascade model |Tl| ], the complexity of 
these calculations makes it difficult to obtain a lucid understanding of the dependence of the different time scales on 
various parameters and model assumptions. It is the goal of our present investigation to derive this insight. 

A second, equally important, motivation for our study was the desire to obtain a better physical understanding 
of the infrared cut-offs required in the parton cascade model. The two cut-off parameters employed in ref. [1], 
the minimal transverse momentum transfer po in binary parton interactions and the infrared cutoff z m i n in parton 
fragmentation processes, were determined by comparison with cross-sections and particle multiplicities measured 
in nucleon-nucleon interactions at high energies. These values are assumed to reflect the transition between the 
perturbative and nonperturbative regimes of QCD in the normal vacuum. On the other hand, it is well known that 
color screening provides a natural cut-off of long-range interactions in a deconfined QCD plasm a fll2|| , and no artificial 
cut-off parameters are required to obtain finite perturbative cross sections in a dense medium |13|]. 

A previous study p4| has shown that screening effects may be sufficiently strong immediately after the primary 
parton collision events in nuclear collisions at RHIC energies and beyond to provide the infrared cut-off required in 
the treatment of secondary parton interactions. Improved understanding |l5|] of the suppression of fragmentation 
processes in a dense QCD plasma (Landau-Pomeranchuk effect) has also allowed to regard the infrared cut-off z m [ n as 
a medium-dependent effect in relativistic nuclear collisions. We show here how these medium effects can be utilized 
to obtain a parameter-free set of equations, based on perturbative QCD in a dense partonic medium, that describes 
the evolution of quark and gluon distributions towards equilibrium. Arbitrary cut-off parameters enter only into the 
description of the primary semi-hard parton scattering, where one must continue to rely on a comparison with nucleon 
scattering data. After this short initial phase, however, the approach towards an equilibrated quark-gluon plasma will 
be described without need for arbitrary parameters. 

Our paper is structured as follows: Chapter II discusses initial parton production in relativistic nuclear collisions. 
The evolution of these partons into locally isotropic, quasi-thermal momentum distributions is described in Chapter 
III. In Chapter IV we derive a set of rate equations describing the further evolution and chemical equilibration of 
gluon and quark distributions. The influence of the dense medium on the relevant QCD cross sections is discussed 
in Chapter V, and a closed set of rate equations including these medium modifications is derived. We explore the 
solution of these equations in Chapter VI. Some consequences for quark-gluon plasma signatures are briefly discussed 
in the Conclusions. 



II. MINUET PRODUCTION 



In high energy nucleon-nucleon collisions the production of minijets with pt about a few GeV becomes increasingly 
important at colliding energies beyond the CERN ISR energy range jl6| . One would certainly expect that there could 
be a fairly large number of minijets produced in ultrarelativistic heavy ion collisons. It was estimated by Kajantie 
et al. fl7| l that minijets could contribute up to half of the total transverse energy produced in heavy ion collisions at 
RHIC energy. Unlike the soft processes which dominate the reaction at lower energies, those minijets carry relatively 
large transverse momenta and they quickly become incoherent from the rest of nuclear matter in the fragmentation 
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region. If a sufficiently high density is reached and the equilibration proceeds rapidly, those minijets would eventually 
lead to a thermalized and locally deconfined quark gluon plasma. 

To estimate the initial parton density, we calculate the minijet production using the HIJING Monte Carlo model 
Irsf . In this model, perturbative QCD, implemented along the lines of Pythia (T9), is combined together with low 
Pt phenomenology to describe multiple minijet production in each binary nuclcon-nucleon collision. One important 
quantity in the model is the inclusive jet cross section. Given the parton structure function, f a (x,p T ) and the 
perturbative parton-parton cross section do~ a b, the differential jet cross section in nucleon-nucleon collisions can be 
calculated as M 



da 



jet 



dp T dyidy 2 



= K^2x 1 f a (x 1 ,p^)x 2 fb(x2,PT) 



a , b 



dOgb 

di 



(1) 



where the phenomenological factor K ss 2 accounts for higher order corrections and x\ , xi are the Feynman vari- 
ables denoting the longitudinal momentum fraction carried by a parton. Since this cross section diverges and the 
perturbative approach fails at small pt, we introduce an infrared cut-off po to calculate the total inclusive jet cross 
section a; e ^(po). Introducing a so ft for the soft interactions below the cut-off, we obtain the inelastic cross section for 



nucleon-nucleon collisions in the eikonal approximation 21 



<r m = I d b 



1-e 



T soft +ff jet> T ~W 



(2) 



where Tjv (b) is the partonic overlap function between two nucleons at impact parameter b. In order to account for the 
shadowing of the parton density within a nucleus, we also introduced effective parton structure functions f a /A(x,p T ) 
which include medium corrections p"8[ |. 

To fix the two correlated parameters po and c g0 ^, we have assumed constant values of po — 2 GeV/c and cr g0 ^ 
= 57 mb. The resultant total, elastic and inelastic cross sections agree very well with the experiments from ISR 
to Tevatron and cosmic-ray energies p2| . We want to emphasize here that the infrared cut-off po in the primary 
parton-parton scatterings is a phenomenological parameter and is constrained by the experimental values of the total 
nucleon-nucleon cross section. (As we will discuss below, this infrared singularity of parton cross sections in a dense 
partonic medium can be regulated naturally by the color screening mass.) 

Due to the high gluon density at small x and the large gluon-gluon cross section, the initially produced partons are 
predominantly gluons. It is very important for the following discussion of the emergence of momentum isotropy that 
the pt distributions of the initially producedgluons are almost exponential. The results of the HIJING calculation 
at RHIC and LHC energies are given in ref. H|, where the rapidity width of the central plateau Y, the total number 
of the produced gluons Nq, and their average transverse momentum {kx) are listed. Before we estimate the initial 
parton number and energy densities, we next have to address the problem of thermalization. 



III. THERMALIZATION 



We can apply rate equations and particle distribution functions to describe the system only, when approximate local 
isotropy in momentum space is achieved. At the onset of expansion, scattered partons with very different rapidities 
are confined to a highly compressed slab. The width A ss 2/po = 0.2 fm of this slab is determined by the virtuality of 
the scattering partons, which in turn is connected with the transverse momentum cut-off po for perturbative parton 
scattering. 

At this point we adopt a realistic ansatz for the phase space distribution of the scattered gluonic partons. Following 
ref. we define the momentum distribution function as: 

where kx and y are the transverse momentum and the rapidity of the gluons respectively, and d& — 16 is the color-spin 
degeneracy factor for gluons. Parametrizing the HIJING results g(kx,y) can be represented by 

9{kT,y) = ^[&(y + Y)-Q{y-Y)] g(k T ), (4) 
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where the function g(kx) is explicitly given in ref. H] for Au + Au collisions at RHIC and LHC energies. For the 
spatial part of the distribution function, we assume that the partons form a homogeneous cylinder of length A and 
radius R. The phase space distribution function then factorizes: 



with D (x) given by 



The overall normalization 



F(k,x) = f(k)D(x), (5) 



Nq = J d 3 k J d 3 x F(k,x) = J d 3 k f(k) (7) 

yields gluon rapidity densities N G /2Y = 110 for RHIC and N G /2Y = 810 for LHC. 

By local isotropy in momentum space we mean that a parton located in the central region of the collision is 
surrounded by partons with different but isotropically distributed momenta. This parton does not interact with all 
the partons in the plasma, but only with those within a distance approximately equal to the parton mean free path. 
We can imagine a little box in the central region the size of the mean free path A f . The momentum of each particle in 
the box has components longitudinal and transverse to the beam, kh and kx respectively. We say that local isotropy 
in momentum space is established when the variance <jl of the longitudinal momentum distribution of partons in the 
box equals that of the transverse momentum distribution, <tt- In the case of an exponential distribution we call the 
isotropic parton system (approximately) thermalized. The variance <7t is defined as 

4 = \ (k 2 T ) = ^J d3k k rf ( fc ) ■ (8) 

At RHIC energy we find <tt = 1-07 GeV/c using the phase space distribution of ref. [[uj, the value at LHC energy is 
1.76 GeV/c. For the mean free path Ay of this initial phase we take pij : 

3 

A" 1 = -a s a T , (9) 

corresponding to Af — 0.4 fm (RHIC) and Af = 0.25 fm (LHC), assuming a s = 0.3. 

As the system expands, mainly in the longitudinal direction, partons with longitudinal momentum will leave the 
central region and <jl will decrease as a function of time. Assuming that the system evolves approximately by free 
streaming, the longitudinal position of a parton depends on its rapidity and time 

z(t) = 2(0) + tanh(y)i. 

Here t denotes time in the cm. system and z(0) is the initial position of the parton along the beam axis (|z(0)| < 
The phase space distribution F(k, x) depends accordingly on time, 

F(k, x; t) = F(k, x - tanh(y)tz ; 0), (10) 

where z is a unit vector along the ^-direction. Now we can determine and ut as functions of time: 

2 (f s Jd 3 xJd 3 kF(k,x;t)kl 

L{, ~ Jd 3 xJd^kF(k,x;t) ' 1 ' 

where the x integral is over Vf, the volume of a tube of longitudinal length Af, centered at z = 0. We find that <JL(t) 
equals <7t for tj so = 0.31 fm/c (RHIC) and t^ SQ = 0.23 fm/c (LHC). For t > ij g0 , there is local isotropy in momentum 
space and we will describe the system by thermal rate equations from then on. 

IV. PARTON CHEMISTRY 

As discussed in the previous section, the kinematic separation of partons with different rapidity establishes conditions 
required for the validity of continuum dynamics after a short time of the order of 0.3 fm/c. At this time the momentum 
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space distribution of partons is roughly isotropic locally and approximately exponential. Since we are here primarily 
interested in the chemical equilibration of the parton gas, we shall assume that the parton distributions can be 
approximated by thermal phase space distributions with non-equilibrium fugacities A^: 



f(k; T, A,) = A?e-"^ fc = A, (e^ fe ± A,) 1 , (12) 

71=1 

where [3 is the inverse temperature and u M is the four-velocity of the local comoving reference frame. The expression 
(O) is known as Jiittner distribution. When the parton fugacities \ are much less than unity as may happen during 



the early evolution of the parton system, we can neglect the quantum corrections in (12) and write the momentum 
distributions in the factorized form 

f(k;T,\ i )*X i e-> 3u - k . (13) 

However, this introduces errors of the order of 40 percent, when the parton phase space distributions approach 
chemical equilibrium. As an example, we show the dependence of the screening length fj.^ 1 for static color-electric 
fields related to the square of the Debye mass Mm 



^2 

,2 



6q 2 f°° 

4r / kf(k)dk, 
^ Jo 



fijj = -V / kf(k)dk, (14) 



on the gluon fugacity X g in Figure 1. The solid line is the exact perturbative result for the Jiittner distribution 
(|l2l), the dotted line shows the Boltzmann approximation (|l3|), and the dashed line shows the result obtained for the 
factorized Bose distribution 

f(k;T,\i) =A i (e" u - fe ±l) _1 . (15) 

The dimensionless quantity (ij-,/ X g g 2 T 2 drops from unity at X g — 1 to 6/ir 2 at 
A g = 0. As can be seen the Boltzmann distribution function ([l3]) is a good approximation to the gluon phase 
space density at small A g w values. For A g w 1 values, however, the approximation ( [l5|) applies, which we will 
adopt in most of the following calculations. Similar, although smaller, violations of factorization occur in the gluon 
and quark densities and in the various reaction rates. 

In general, chemical rections among partons can be quite complicated because of the possibility of initial and final- 
state gluon radiation. Multiple radiative emission processes have been investigated in the framework of numerical 
simulations of parton cascade s El BI. However, since radiation processes are considerably suppressed in a dense parton 
medium due to rescattering pa], we shall here consider only processes where a single additional gluon is radiated, 
such as gg — > ggg. In order to permit approach to chemical equilibrium, the reverse process, i.e. gluon absorption, 
has to be included as well, which is easily achieved making use of detailed balance. Closer inspection shows that gluon 
radiation is dominated by the process gg — ► ggg, because radiative processes involving quarks have substantially 
smaller cross sections in perturbative QCD, and quarks are considerably less abundant than gluons in the initial 
phase of the chemical evolution of the parton gas. Here we are interested in understanding the basic mechanisms 
underlying the formation of a chemically equilibrated quark-gluon plasma, and the essential time-scales. We hence 
restrict our considerations to the dominant reaction mechanisms for the equilibration of each parton flavour. These 
are just four processes [ pq ]: 

99 *-> 999, 99 <-> QQ- (16) 

Other scattering processes ensure the maintenance of thermal equilibrium (gg «-> gg, gq <-> gq, etc.) or yield corrections 
to the dominant reaction rates (gq «-> qgg, etc.). 

Restricting to the reactions (Hq) and assuming that elastic parton scattering is sufficiently rapid to maintain local 
thermal equilibrium, the evolution of the parton densities is governed by the master equations p5fl : 

d^ngU^) = n g (i? 2 ^3 - #3^2) - (n g R g -> q - n q R q ^, g ) (17) 

d^UgU^) = d^UqU^) = n g R g ^ q - TlqRq^g, (18) 

where i?2^3 and i?3^2 denote the rates for the process gg — > ggg and its reverse, and R g ^ q and R q ^ g those for the 
process gg — > qq and its reverse, respectively. The temperature evolves according to the hydrodynamic equation 

9„(e/) + Pa/=0, (19) 
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where viscosity effects have been neglected |26j . 

The rates depend on the temperature and the fugacities in a rather complicated way. In order to make analytical 
progress, we here assume factorization of the dependence of the densities on the fugacities, i.e. eq. (|l5|), and postpone 
a full numerical evaluation of the rate equations. The rates can then be factorized in the following simple way: 

n g (i? 2 ^ 3 - R 3 ^ 2 ) = ±a 3 n 2 g ( 1 - , (20) 



n 



9 I IL q ll! q lb n \ 

UgRg^g - U q R q ^ g = O 2 U g 1 - (21) 

where 03 and a 2 are thermally averaged cross sections: 

0-3 = (v(gg ggg)}, cr 2 = (<r(gg -» qq)) (22) 
and fii denote the densities for Aj = 1. The equations ( 17|l8| ) can then be rewritten as 

d^XghgV?) « n g R 3 X g {l - \ g ) - 2n g R 2 X g (l - , (23) 

d^XghgU") = d^XghgU") « h g R 2 X g ^1 - , (24) 

where the density weighted reaction rates R 3 and R 2 are defined as 

#3 = |<^3« ffJ R2 = h^2n g . (25) 



In order to obtain simple solutions to eqs. (|2^ 
longitudinal, yielding Bjorken's scaling solution 



24j) we will assume that the expansion of the parton fireball is purely 
of the hydrodynamic equation (|l9|): 

^ + ^ = 0, (26) 
ar t 



where r is the proper time. This assumption is expected to be very well satisfied during the early expansion phase 
of the fireball, especially at proper time r -C Ra, where Ra is the transverse radius of the fireball. Further assuming 
the ultrarelativistic equation of state 

e = 3P= [a 2 X g + b 2 (X q + Xg)] T 4 , (27) 

with a 2 — 87r 2 /15, b 2 — 77r 2 iVy/40, where Nf is the number of dynamical quark flavors, we obtain from ( p^ ) a relation 
between the proper-time dependence of the temperature and the fugacities. Denoting the derivative with respect to 
t by an overdot, we find 

^±M +4 f * 0; (28) 
X g + b(X q + Xg) T 3r 

where 

b = b 2 /a 2 = ^l. (29) 

Equation ( p8| ) is easily integrated, resulting in 

[X g + b(X q + Xg)] 3/4 T 3 T = const. (30) 

For a fully equilibrated quark-gluon plasma (X g = X q = Xq = 1) this corresponds to the Bjorken solution T(r) = 
ro(ro/T)V3. 

In the same approximation, i.e. neglecting interactions and using the factorized distributions, the gluon and quark 
equilibrium densities are related to the temperature parameter T by: 



n g = ^C(3)T 3 = ai T 3 , (31) 



16 

^((3)N f T 3 = hT 3 . (32) 
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We can also rewrite eqs. (19, 23| , |24| ) as equations determining the time dependence of the parameters T, X g and A g . 
In the longitudinal scaling expansion we find: 



f/77 71 

After division by n„ and n„, respectively, the rate equations take the form: 



3^ + i = # 3 (1 - Xg) ~ 2R 2 



1 _ X q Xq 
*2 



T 



0-1 f X 



x,-, 



with the additional constraint 



h \X 



(\ ~ A,)T 3 



const. 



(33) 

(34) 
(35) 

(36) 



expressing baryon number conservation. Being interested in parton thermalization at very high collision energies, we 
shall assume baryon symmetric matter, i.e. X q = Xq, which solves @) trivially. Equations (pO|,pM 35) determine the 
evolution of T(r), A s (r), and A 9 (t) towards chemical equilibrium, once the reaction rates R2 and R3 are known. We 
now turn to these. 



V. THE EQUILIBRATION RATES 

The cross sections 02 and 03 contain infrared singularities if calculated in naive perturbation theory. Recently, 
Braaten and Pisarski ]28[ ] proposed an effective perturbation theory, based on a resummation of subsets of diagrams 
(hard thermal loops), which takes screening effects into account and avoids inconsistencies of the naive perturbation 
theory at finite temperature. This method results in using effective propagators and vertices, which show a complicated 
momentum dependence. It has been used, e.g., to calculate the stopping power p9| ] and the viscosity ]3C[ | of the quark- 
gluon plasma, which also follow both from elastic parton collisions. In order to calculate 02 and 173 we will adopt a 
simplified version of this idea by introducing momentum independent screening masses into the propagators, whenever 
infrared divergences arise otherwise. In the case of the stopping power and the viscosity this approximation provides 
quantitatively good results. 



A. Gluon equilibration 

The rate R3 for the process gg — > ggg depends on the triple differential radiative gluon-gluon cross section, which 
can be written as [Bf| 



da 



2^3 



da 



2^2 



d 2 q±dyd 2 k± d 2 q± 



Caol s 



(37) 



where Ca = 3 is the Casimir operator of the adjoint representation of SU(3). Here denotes the transverse 
momentum and y the longitudinal rapidity of the radiated gluon, and denotes the momentum transfer in the 
elastic collision. The in-medium cross section for elastic scattering of two gluons 64] 



da 2 



d 2 qi 



C a 



2a 2 



(38) 



with C gg = |, is screened by the Debye mass. 

In the presence of a dense medium the emission of radiation is suppressed, if the gluons scatter again before the 
emission is completed (Landau-Pomeranchuk effect), leading to the condition p3[ 



k±Af > 2 coshy, 



(39) 



where Ay is the mean free path of a gluon. The contribution from soft radiation is strongly suppressed by this effect. 
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To obtain the gluon production rate R3 we must integrate the differential cross section ( |37| ) over momentum transfer 
q ± and the phase space of the radiated gluon. The integrand is singular at k± = q ± . We shall deal with this problem 
in two different ways. First, we shall assume that k± ■ q ± = in the integrand, i.e. that the radiation is mostly 
perpendicular to the reaction plane. Later we will average over the angle between k± and introducing an infrared 
cut-off to avoid the singularity. In the first approach, we replace (fej_ — q±) 2 by (kj_ + q±) in the integrand, effectively 
suppressing collinear bremsstrahlung. We then arrive at the following modified differential cross section: 

//n-2-»3 ICaCI r\ 3 n 2 



where the step functions account for the Landau-Pomeranchuk effect (39) and for energy conservation, (s) = 1ST 
the average squared center-of-mass energy of two gluons in the thermal gas. The integrated elastic gluon-gluon cross 
section in the medium is from eq. 



a 2 ^ 2 = 2nC gg ^-, (41) 

yielding a fugacity independent mean free path 

A] 1 = (J 2 ^ 2 n g = fai a s T. (42) 

Using these values we evaluate the chemical gluon equilibration rate R3 = ^n g as, as defined in eq. (gij), numerically. 
This rate scales with the temperature linearly but is a complicated function of the gluon fugacity. The solid line in 
Figure 2 shows the scaled rate R3/T versus X g for a coupling constant a s = 0.3. The dotted line corresponds to the 
analytical fit 

R 3 = 2.1a 2 s T(2\ g -\ 2 g ) 1/2 , (43) 

which will be used in solving the time dependent rate equations discussed in the next section. 

We have also studied another solution to the problem posed by the collinear singularity in eq. ( |37| ) , namely, to add 
a screening mass to (k± — q±) 2 - The integral over the relative angle (p between k± and q ± can then be performed 
analytically, yielding 

± J d<p[(k ± - q ± f + M 2 ]- 1 = Pi + q\ + M 2 ) - 4fci,i]- 1/2 . (44) 

The rate R3 obtained in this way for fi = [ip is shown in Figure 2 by the dashed line. It is somewhat larger than the 
rate obtained by setting cos^ = 0, so that (E3I) can be considered as a conservative estimate. 



B. Quark equilibration 

The total cross section for the process gg — > qq is dominated by the Compton diagrams involving exchange of a 
virtual light quark in the t- or u- channel. For massless quarks the differential cross section |32]] 



dt 



Twq 

a 2 



"1 (u 




3 


3ut 


[e [I" 


-=) 


8 


h 4^_ 



(45) 



diverges as u, t — > 0, hence the medium-induced effective quark mass plays a crucial role. Unfortunately, a complete 
calculation of the cross section for the process gg — » qq within the framework of thermal field theory with hard thermal 
loop resummation |^8| has not yet been done. But since the divergence of the cross section is only logarithmic in 
this case, it may be a sufficiently good estimate to simply substitute the effective thermal quark mass as cutoff in the 
di verg ent integral over momentum transfer. The thermal quark mass, using Bose and Fermi equilibrium distributions 
is 



(46) 



yielding the total cross section (for s » M 2 ): 
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s 7 



3s 

Integrating over thermal gluon distributions and inserting the average thermal (s) = 1ST 2 in the logarithm, we have 

For a s = 0.3 in the logarithm and neglecting X q , we obtain the light quark production rate 

i? 2 = \o2n g « 0.064/V, c^A ff T ^ln . (49) 



VI. RESULTS 



We have solved the rate equations (34 3^) and the energy conservation equation ( p8| ) simultaneously by numerical 
integration using a fourth order Runge-Kutta method. The initial conditions for these rate equations are the number 
density no = n(ij g0 ) and the transverse energy density et — £ T{t[ S0 ) of gluonic partons, where 

1 N G 

no = ^m- W ^ = n (k T ). (50) 

The total initial gluon energy density is £o = Aex/it. We obtain no = 2.4 fm -3 and £o = 3.5 GeV/fm 3 at RHIC 
energy, and no = 23 fm~ 3 and e = 52 GeV/fm 3 at LHC energy. Quarks contribute a smaller amount to the initial 
parton energy density, because the quark-producing cross sections are smaller in perturbative QCD than those for 
gluon production p^| . The total quark contribution to the energy density we estimate as 30 percent, yielding total 
initial energy densities of 5 GeV/fm 3 at RHIC and 70 GeV/fm 3 at LHC. 

The temperature and fugacities at the beginning of hydrodynamical evolution (t = ij so ) are obtained assuming a 
thermal, but not chemically equilibrated Bose distribution for gluons, i.e. 

n = A« ai T 3 , £ = A^ 2 r 4 (51) 

where a\ m 1.95, a<i » 5.26. The quark fugacity is taken as A° = A°/5, corresponding to a ratio 3:1 of the initial gluon 
number to the combined number of quarks and antiquarks. The Table shows these relevant quantities at the moment 
t[ so , for Au + Au collisions at RHIC and LHC energies. 

The evolution of temperature and the fugacities are shown Figures 3. We find that the parton gas cools considerably 
faster than predicted by Bjorken's scaling solution (T 3 r = const.), because the production of additional partons 
aproaching the chemical equilibrium state consumes an appreciable amount of energy. This can be also deduced from 
eq. (pOf): while the fugacities increase, the product T 3 r must decrease. 

The accelerated cooling on the other hand impedes the chemical equilibration process, which is more apparent at 
RHIC (Fig. 3a) than at LHC energies (Fig. 3b). In order to see where the perturbative description of the parton 
plasma is applicable we investigate the time evolution of the total energy density, eq. (|2^). The solid lines in Fig. 4 
correspond to the initial conditions shown in Table I, while the dashed line corresponds to a more optimistic estimate 
of the initial gluon production at RHIC. One realizes that the perturbative parton plasma has a lifetime of 1 — 2 fm/c 
at RHIC, while at LHC the plasma may exist in a deconfined phase as long as 5 — 6 fm/c. The evolution of the energy 
density is remarkably insensitive to the magnitude of the gluon equilibration rate R3. A change in i?3 by a factor 2 
is hardly noticeable in the function e(r). 

The initial conditions obtained in section III depend sensitively on the low-momentum cut-off po used to regularize 
the mini-jet cross sections. This parameter has usually been determined from fits to the inelastic nuclcon-nucleon 
cross section. Ambiguities of this method allow for a range of 1.4 — 2 GeV for the low-momentum cut-off at RHIC 
energy. Using the smaller value [ p5[ one obtains an initial gluon fugacity X g m 0.18, which is three times larger than 
the value discussed above. Using this value we find a larger initial energy density and hence a longer lifetime of 2.5 
fm/c for the plasma at RHIC (see the dashed curve in Fig. 4). A higher value of X g at t^ SQ can also result from 
rescattering of gluons during the "free-streaming" period between the end of initial semihard scattering and ij so - 
Indeed, the parton cascade simulations [9 indicate a rapid increase of gluon density during this pre-thermal phase. 
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We conclude that the validity of this scenario of quark-gluon plasma production by perturbative parton scattering is 
uncertain at RHIC energies, but very likely in the LHC energy domain. 

From our investigation emerges the following scenario of a nuclear collision at collider energies: Within 0.2 — 0.3 
fm/c a dense parton gas is produced at central rapidities, which can be described as a locally thermalized, but not 
chemically equilibrated quark-gluon pasma. In the fragmentation region, where the parton gas created by minijets is 
not sufficiently dense to screen color fields, strings and color ropes are formed |p6| , p7| , which decay by non-perturbative 
QCD processes on a time scale of about 1 fm/c. These strings extend between the wounded nuclei and the surface of 
the parton plasma, penetrating up to a screening length determined by the actual Debye mass. 

As in the Bjorken scenario, the scaling expansion of the system can best be described by using the proper time r 
and space-time rapidity 

ry = Iln^ (52) 

t — x 

variables. This scenario is schematically shown in Fig. 5, where the physically different space-time regions of the 
pre-thermal creation phase, the parton plasma, and the string-dominated phase are indicated. 

The main difference to the conventional Bjorken scenario p7| is the far shorter formation time in the central rapidity 
region, corresponding to a much higher initial temperature. The Bjorken scenario, however, most probably applies to 
the fragmentation region, where the non-perturbative color ropes may also decay into a quark-gluon plasma. At the 
interface between the two regions a reheating of the parton plasma may occur by dissipation of field energy stored in 
the strings. The dissipation of field energy in a quark-gluon plasma was recently studied by Eskola and Gyulassy [ g8| . 
At RHIC energies, where the approach discussed here leads to a rather low initial density of partons, color flux tubes 
may still permeate the whole central region. This would result in multiple interactions between partons from minijet 
cascades and the soft component modeled by color strings or ropes. In addition, the temperature falls so rapidly 
that it is difficult to maintain thermal equilibrium by elastic two-body collisions among partons. At RHIC energies 
the characteristic collision rate, given by the inverse mean free path A7 , is comparable to the cooling rate T/T. As 
discussed by the authors of ref. [38] , such a state would exhibit large viscosity effects which can substantially slow 
down the cooling rate of the QCD plasma. 



VII. CONCLUSIONS 



Our results support the picture emerging from numerical simulations of parton cascades and simple estimates that 
the highly excited parton plasma created in ultrarelativistic nuclear collisions is initially mainly a gluon plasma. Our 
results have important implications for several experimental signals associated with quark-gluon plasma formation, 
as has been pointed out before [8]- pd|j39| , ^fH : Rapid gluon thermalization at a high initial temperature leads to a 
substantial thermal contribution to the total yield of charmed quarks. This increase comes on top of a substantial 
suppression of primary charm production by nuclear shadowing of gluon distributions in the central rapidity region. 
Hence the unambiguous observation of this charm enhancement will require the comparative study of p + p, p + A, 
and A + A collisions at the same center-of-mass energy. 

The lack of chemical equilibration in the quark and gluon densities, on the other hand, causes a severe depletion 
of the number of emitted lepton-pairs compared with the naive thermal estimate. However, the shift towards higher 
invariant masses of the steeply falling lepton-pair spectrum due to the larger initial temperature in the plasma can 
conceivably offset this suppression J39| , p0| | . 

Reliable quantitative predictions for these experimental signatures require complete microscopic simulations of the 
parton cascade which include the medium modification of cross-section and effective masses discussed in this paper. As 
we have shown, these medium effects slow the chemical equilibration of the parton plasma down. How can the medium 
dependence of the parton cross sections be included in microscopic simulations? Since these calculations trace the full 
space-time evolution of parton densities, it would be straightforward to calculate the color screening length, as well 
as the effective medium-induced masses of collective gluon and quark modes (plasmons and plasminos) as function of 
time. These quantities could then be used as input parameters for the determination of medium-dependent scattering 
cross sections and branching rates. 

An unresolved issue remains in the medium-dependence of the lower cut-off po used to regularize perturbative QCD 
cross sections, which was fixed at po — 2 GeV/c in HIJING and at a similar, but energy-dependent value in the parton 
cascade. One possible way of approach would be to identify an upper limit to the (medium-dependent) running QCD 
coupling constant, above which perturbative cross sections cannot be used. This limit could be chosen such 

that it equals a s (po) in vacuum; but it would correspond to a different momentum cut-off, or even no cut-off at all, 
at high parton densities p| . 
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RHIC 


LHC 


*iso ( fm / c ) 


0.31 


0.23 


eo (GeV/rrn J 


3.5 


52 


no (fm~ 3 ) 


2.4 


23 


etot (GcV/fm 3 ) 


5 


70 


A / ( fm ) 


0.40 


0.24 


(fc T > (GeV) 


1.17 


1.76 


T (GeV) 


0.55 


0.83 


x° 

A 9 


0.06 


0.16 


A, 


0.012 


0.03 



TABLE I. Values of the relevant parameters characterizing the parton plasma at the moment tj so , when local isotropy of 
the momentum distribution is first reached. 
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FIG. 1. The square of the thermal screening mass no of color electric fields scaled by the equilibrium value g 2 T 2 is depicted 
as function of the gluon fugacity A 9 . The solid curve shows the correct result for the Jiittner distribution ([12]), the dashed 
line corresponds to the factorized Bose distribution (|l5|), while the dotted line indicates the Boltzmann approximation ( |l3| ) 
applying in the limit of a very dilute gluon plasma. 

FIG. 2. The scaled gluon production rate R3/T is shown as function of the gluon fugacity A 9 for a s — 0.3. The solid line 
shows the numerical result, using the approximation (42), while the dotted line shows the analytical fit (45). The dashed line 
represents the result obtained with the screening approximation (46) to the collinear singularity in the bremsstrahlung cross 
section. 

FIG. 3. Time evolution of the temperature T and the fugacities A 9 and X q of gluons and quarks in the parton plasma created 
in Au + Au collisions at (a) the RHIC cm. energy of 100 GeV/nucleon, (b) the LHC cm. energy of 3000 GeV/nucleon. The 
initial values for T, X g and \ q are determined from simulations with the HIJING code and are listed in the Table. 

FIG. 4. Time dependence of the energy density of the parton plasma at RHIC and LHC energies. The solid lines show the 
results obtained with the initial values taken from the Table. The dashed line shows the effect of a three times higher initial 
gluon density at pRHIC energy, corresponding to QCD cut-off po ~ 1.5 GeV instead of po = 2 GeV. Whereas a well-developed 
parton plasma is predicted to be formed at the LHC, the prediction for RHIC is uncertain. 

FIG. 5. Schematic space-time diagram of the evolution of QCD matter in the initial phase of an ultrarelativistic heavy 
ion collision. The central rapidity region is dominated by a thermalized, but not chemically equilibrated parton plasma after 
t iso w 0.2 - 0.3 fm/c. 
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